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We present an experimental study of density and order fluctuations in the vicinity of the solid- 
liquid-like transition that occurs in a vibrated quasi-two-dimensional granular system. The two- 
dimensional projected static and dynamic correlation functions are studied. We show that den- 
sity fluctuations, characterized through the structure factor, increase in size and intensity as the 
transition is approached, but they do not change significantly at the transition itself. The dense, 
metastable clusters, which present square symmetry, also increase their local order in the vicinity of 
the transition. This is characterized through the bond-orientational order parameter Q4,, which in 
Fourier space obeys an Ornstein-Zernike behavior. Depending on filling density and vertical height, 
the transition can be of first or second order type. In the latter case, the associated correlation 
length ^4, relaxation time r4, zero k limit of Q4 fiuctuations (static susceptibility), the pair cor- 
relation function of Q4,, and the amplitude of the order parameter obey critical power laws, with 
saturations due to finite size effects. Their respective critical exponents are = 1, u\\ = 2, 'y = 1, 
rj = 0.67, and P — 1/2, whereas the dynamical critical exponent z — i'^^/i'± — 2. These results are 
consistent with model C of dynamical critical phenomena, valid for a non-conserved critical order 
parameter (bond-orientation order) coupled to a conserved field (density). 

PACS numbers: 64.60.Ht, 45.70.-n, 05.40.-a, 64.70.qj 



Introduction. A dry granular system is an a-thermal 
collection of macroscopic particles that interact mainly 
through dissipative, hard core-like, collisions. They be- 
have as solids, liquids or gases depending on the nature 
of the forces that act upon them and the energy injec- 
tion rate [if. These systems present phase transitions 
and coexistence. Simple examples are: a thin layer of vi- 
brated sand, which for small driving amplitudes remains 
solid but can be completely fluidized for larger acceler- 
ations; an avalanche driven by gravity were a thin layer 
of grains flows above an almost solid pile. Granular sys- 
tems are therefore excellent candidates for studying non- 
equilibrium phase transitions 0j 0] • In vibrated thin lay- 
ers energy is transferred from the top and bottom lids to 
the vertical motion of the grains, which later transfer the 
energy to the horizontal motion at collisions, that arc 
also dissipative. The sequence breaks detailed balance 
keeping the system out of equilibrium. 

Recently, several granular systems that undergo in- 



teresting phase transitions have been reported |4|-|15 



One particular system is a vibrated fluidized granular 
monolayer composed of N hard spheres of diameter d 
confined in a shallow cell of height Lz < 2d (typi- 
cally Lz sa 1.7d— 1.9d). Under proper conditions, solid 
and liquid phases can coexist at mechanical equilibrium 
(3.[l0l U, 14|. The solid clusters can present different or- 
der symmetries, like square or hexagonal, depending on 
forcing, geometrical and particle parameters. It has been 
reported that for « 1.7d — 1.8d and for a large range 
of filling densities, the most compact structure in quasi- 
2D is made of two layers of square symmetry. The more 
compact hexagonal structure formed by two layers needs 
a larger vertical gap or larger densities [ll|. The criti- 



cal amplitude above which there is coexistence decreases 
with increasing density. 

Many of the previous works on granular phase transi- 
tions focus on the similarities or comparisons of such non- 
eq uilibrium sy stems with equilibrium phase transitions 
For example, the equilibrium KTHNY 



lOhUJ, 15 



theory has proved useful in the two dimensional melt- 
ing of granular monolayers [l2j . Here we focus on a 
dynamical critical phase transition in a non-equilibrium 
quasi-2D granular system. We present an experimental 
study of the solid-liquid phase transition in a vibrated flu- 
idized granular monolayer. The solid phase consists on 
two square interlaced layers, stabilized by the collisions 
with the top and bottom walls and the confining pressure 
exerted by the liquid phase ll|, . Wc focus on den- 



sity and bond-orientation order fluctuations in the vicin- 
ity of the transition. We show that the transition can 
be continuous or abrupt depending on the cell's height 
and fllling density. Density fluctuations show a crossover 
behavior at the transition, whereas order shows strong 
fluctuations. In the continuous case several magnitudes 
show critical-like behavior being possible to measure flve 
independent critical exponents. These results are consis- 
tent with model C of dynamical critical phenomena [lij , 
valid for a non-conscrvcd critical order parameter (bond- 
orientation order) coupled to a conserved field (density) . 

Experimental setup and procedures. The experimen- 
tal device is similar to the one described in [17|. The 
granular system is composed by ~ 10^ stainless steel 
spherical particles, of diameter d = 1 mm. The quasi- 
two-dimensional box has lateral dimensions = Ly = 
L ~ lOOd. In order to study the different possible 
transitions, first or second order type, two configura- 
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FIG. 1. (color online), (a) S{k) in the large wavelength limit for four different accelerations, F = 2.13, 2.74, 3.18 and 4.79 (C2). 
Open symbols correspond to raw S{k) data, whereas solid symbols correspond to averages using windows k £ [nfcmin, (?^+l)fcmin] 
for integer n ^ 1, where femin = tt/L. Error bars correspond to standard deviations. The inset shows S{k) for a larger range of 
k for F — 3.64. (b) Pre-peak maximum Smax = S{k — k*), which occurs at fc = fc*, and (c) associated length scale ^ — Tr/k* as 
functions of F for CI (open symbols) and C2 (closed symbols). 



tions are used: = 1.83d ± 0.02rf and N = 9878, 
or L, = IMd ± Omd and TV = 11504. These are 
labeled configurations 1 (CI) and 2 (C2) respectively. 
The surface coverage is defined by the filling fraction 
(j) = Nnd^/AL^ {(l)ci = 0.776 and (t>c2 = 0.904). 

The whole setup is forced sinusoidally with an 
electromechanical shaker, with displacement z{t) ~ 
AsiTi{ujt). Top view images are obtained with a cam- 
era at 10 fps. Particle positions are determined at sub- 
pixel accuracy. Results have been obtained by fixing the 
particle number N , cell height and driving frequency 
/ = a;/27r = 1/T = 80 Hz. The dimcnsionless accelera- 
tion r = Au]"^ / g is varied in the range 1 — 6. 

Static structure Junction. Particle positions rj (t) in 
the plane (x, y) are determined for each time t. Exper- 
imentally, there is no access to the z coordinate. Thus, 
the 2D microscopic density field Fourier components are 

N 



p{k,t) 



d\e^^^-^p{r,t) = Y.^ 



:ik-rj(t) 



(1) 



The static structure factor S{k) measures the intensity 
of density fiuctuations in Fourier space: 



S{k) 



(|p(fc,t)-(p(fc,t))|- 
N 



(2) 



where ( ) denotes time averaging. In general {p{k)) ^ 0, 
due to inhomogeneities induced by boundary conditions. 
The wave vectors are computed from k = TT{nxi + nyj)/ L, 
where n^, Uy S N. 

In the liquid phase (F < Fc) we have verified that 
the system is isotropic, iS'(fc) = S'(fc), where = k. 
In the phase separated regime this quantity is not well 
defined as density fluctuations should in principle behave 
differently in each phase. However, as it is not easy to 
determine Fc precisely, we use this quantity as a possible 
relevant measurement in the vicinity of the solid-liquid 
phase separation, even above the critical amplitude. 



Figure [T^ presents S{k) obtained for C2 (the qualita- 
tive features are the same for both configurations and 
their differences will pointed out explicitly later). The 
main figure presents the long- wavelength range, kd ^ 1, 
for four F below Fc. The inset presents S{k) for a larger 
range of fc. It has the usual form expected for liquids 
with short range order, but with a pre-peak located in 
the range kd = 0.1 — 0.3. The associated density fluctu- 
ations arc indeed visible by simple visual inspection [l^l . 

The pre-peak is characterized by its maximum value 
at k*, S'niax = S{k*), and the associated characteris- 
tic length scale ^ = ir/k*. These quantities are plotted 
in Fig. [T)3-c as functions of F for increasing amplitude 
ramps and for both configurations. They both increase 
as the transition is approached. The difference between 
conflgurations is mainly manifested in the shape of each 
curve, being their final values (near the transition) very 
simflar, S'niax ~ 0.5 — 0.8 and £,/d ~ 20 — 30. By observing 
visually the persistence of the solid clusters we conclude 
that for CI the transition is located at Fc « 2. For C2 
it is more difficult to determine with the same precision 
but it is found to be Fc ^ 5. However, neither S'niax or ^ 
show evident changes at these values. 

Density fluctuations do not show critical behavior, but 
they are needed to create regions of high order. Similar 
density fluctuations have been observed in amorphous 
materials [l^ H^l , which have been consistently related 
to the existence of medium-range-crystalline-order. In 
our case, medium range order will be analyzed with an 
appropriate bond-orientational order parameter, which 
presents critical behavior. 

Bond-orientational order parameter. In the vicinity 
of the transition, fluctuations of high density present the 
same square symmetry as the solid phase. In the quasi- 
2D geometry the solid phase consists of two square inter- 
laced layers instead of the hexagonal layer that is charac- 
teristic of 2D systems [llj . The local order can be charac- 
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terized through a 4-fold bond-oricntational order param- 
eter. This is still valid in quasi-2D geometry because the 
interlaced two-layer square lattices (with unit cell length 
d in each plane) result also in a square lattice when pro- 
jected in 2D, with unit cell length \/2d/2 when the grains 
are close packed. The 4-fold bond-orientational order pa- 
rameter per particle is defined [IJ] 

1 

^ s=l 

where Nj is the number of nearest neighbors of particle 
j and ai is the angle between the neighbor s of particle j 
and the x axis. For a particle in a square lattice, \QW = 1 
and the complex phase measures the square lattice ori- 
entation respect to the x axis. The corresponding global 
average and Fourier components are 

/ 1 ^ \ ^ ^ 
= ( ) , Q4k,t) = ^Qie^'^-^^W. (4) 

The average (|(54|) measures the fraction of particles in 
the ordered phase. This quantity is presented in Fig. [2] 
as function of F for both configurations. Results for in- 
creasing (decreasing) F ramps arc represented by open 
(solid) symbols. Two ramp rates arc also reported: slow 
ramps, for which a quasi-static state has been reached, 
and fast ramps, for which it has not. The difference be- 
tween both configurations is evidenced in the jump of 
about 10% that is measured for (|(34|) at the transition 
for CI. Moreover, the position of this jump depends on 
the F ramp rate: for the slow rate the increasing and de- 
creasing ramp jumps coincide, whereas for faster ramps 
the increasing (decreasing) ramp jump occurs at higher 
(lower) F. We use the slow ramp data to obtain a mea- 
surement of the critical acceleration, F^^ = 2.01 ± 0.03. 
By the contrary, the results obtained for C2 show first 
a linear trend for low F and a clear deviation around 
F « 5.1, with no measurable jump. In fact, the deviation 
from the linear behavior obeys a supercritical-like law. 
For F > 5.2 we have fitted the data with the function 
AQ4 = ilQil) ~Qi = c(F - T^Y, where is the ex- 
trapolation of the linear trend observed for lower F. We 
obtain c = 0.029 ± 0.002, FC2 = 5.12 ± 0.01, and the 
exponent of the order parameter amplitude is /3 = 1/2. 
Within experimental errors, the decreasing ramps also 
coincide with the increasing ramps in this configuration 
(isj . Consequently, the transition for configuration CI is 
abrupt, of first-order type, whereas for C2 it is continu- 
ous, of second-order type. 

Local order can also be analyzed through its fluctu- 
ations in Fourier space by means of the 4-fold bond- 
orientational structure factor 

^ (|g4(fc,t)-(g4(fc,i)>P> 
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FIG. 2. (color online). Average global 4-fold bond- 
orientational order parameter (|(54|) versus F for CI (a) and 
C2 (b). Open (solid) symbols represent data obtained for 
increasing (decreasing) F ramps, with the following rates: 
AF/At f» 0.005 min"^ (A, T) and AF/At « 0.02 min"^ 
(o, •). Continuous lines in (b) correspond to fits of the linear 
trend Q4 = aF + 6 for 2.5 < F < 5, with a = 0.011 ± 0.001 
and h = 0.380 ± 0.002, and a supercritical-like behavior 
(IQ4I) =Q]i+ c(r - T^Y, with /? = 1/2, observed for F > 5. 



For both configurations and for F < Fc, Si{k) shows an 
Ornstein-Zcrnikc-likc behavior in the limit fed ^ I [l8l |. 
Siik) « S'4(0)/[l + (^4fc)2], where ^4 and S'4(0) are the 4- 
fold bond-orientational correlation length and static sus- 
ceptibility respectively. 

For configuration CI, the 4-fold bond-orientational 
susceptibility 6*4(0) and normalized correlation length 
^i/d vary weakly as the transition is approached (Fig. 
[5^). Defining the reduced acceleration e ~ (Fc — F)/Fc, 
we obtain that for 0.005 < e < 0.4, 6*4(0) and ^4/c? vary 
in the ranges 0.25—0.5 and 0.5—0.7 respectively. The fact 
that 6*4(0) < 1 and ^4/d < 1 implies that fluctuations of 
the global 4-fold bond-orientational order parameter are 
weak and that there is practically no order correlation 
below the first-order-type transition. 

For C2 the situation is markedly different. Fig. [3^ 
shows that 6*4(0) and S.4/d vary strongly as the transition 
is approached. In the limit e — > they both saturate, pre- 
sumably due to the system's finite size. For e < 3 x 10~^ 
they saturate to 6*4(0) « 20 and ^4/d ss 10 respectively. 
This figure also demonstrates that both quantities follow 
the critical-like behavior, 

64(0) = a£-'', ^4/^ = (6) 

with the critical exponents 7 = 1 and 1/^ = 1. The crit- 
ical divergence with e makes it necessary to fit the Fc 
separately for each case; for details on the procedure see 
the Supplemental Material [l3| • The adjusted critical ac- 
celerations are F^ = 5.09 ± 0.07 and F^ = 5.24 ±0.08 
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FIG. 3. (color online), (a) 5'4(0) (o, •) and (□, ■) versus 
e for CI (solid symbols) and C2 (open symbols), (b) T4/T ver- 
sus e for C2, for three different low k. Continuous lines show 
critical power law fits, with exponents 7 = i^x = I for ^4(0) 
and ^4, and = 2 for n. The fitted critical accelerations are 
= 5.09±0.07, = 5.24±0.08 and = 5.12±0.07 re- 
spectively. The inset of (b) presents T4 versus ^4 for e > 0.03. 
The continuous line shows a fit r4 ~ (^4)^) with the dynamical 
exponent z — i'ii/i'± = 2. 



respectively. Within experimental errors both critical 
accelerations are very consistent, as well as with the 
value obtained from the supercritical behavior of AQ4 
(rC2 = 5.12 ± 0.01). Notice that now these are obtained 
from fits of measured quantities below the transition, 
whereas before, it was obtained with a fit of the order 
parameter above the transition. 

In the hydrodynamic regime < kd <C 1) and 

in the vicinity of the transition, the bond-orientational 
structure factor is theoretically expected to behave as 
Siik) = Coofc"'^"'''. When T = 5.10 a power-law behav- 
ior is indeed observed in the range 0.1 < kd < 1 with an 
anomalous exponent 77 = 0.67 ± 0.01 [l8| . 

As a final evidence of the observed criticality we now 



turn to the characterization of the relaxation time of the 
metastable solid clusters. The relaxation time is com- 
puted through the two-time bond-orientational correla- 
tion function 



Fi{k,T) 



{5Qi{k,t + T)5Qi{k,ty 
N 



(7) 



where * stands for the complex conjugate and SQi{k, t) ~ 
Q4{k,t) — {Qi{k,t)). Our results show that for low 
wavevectors F4{k,T) sa J4(fc, 0) exp(— r/T4(fc)), from 
which the relaxation time T4(fc) is measured. Here, we 
also obtain a critical-like behavior, which is presented in 
Fig. [Sb. The best fit is obtained for n/T = c£"''ii with 
1^1 1 =2, for which the adjusted critical acceleration is 
rP = 5.12 ± 0.07. The relaxation time also seems to 
saturate for small e, which occurs at smaller e for lower 
k, that is for fluctuations of larger size. The inset of Fig. 
[SJd confirms that ~ with a dynamical exponent 

z = = 2. As usual, there is critical slowing down 

in the dynamics. As a consequence, close to the critical 
point, stationary states are obtained after a long relax- 
ation has taken place. Taken that into account, all T 
ramps for C2 are slow. Also, averages are taken for long 
times. 

Critical dynamics. Five critical exponents have been 
obtained from the analysis of the order parameter. In 
the standard notation of critical phenomena these are: 
/? = 1/2, 7 = 1, ?7 = 0.67, z^_L = 1, and z = 2. In equilib- 
rium, the scaling hypothesis predicts relations among the 
critical exponents. It is worth mentioning that the rela- 
tion 7 = (2 — rj)v± is not satisfied, while a -|- 2/3 -I- 7 = 2 
and v^D = 2 — a {D = 2 \s the spatial dimension) can be 
satisfied simultaneously if a = 0. This exponent, asso- 
ciated in equilibrium to the specific heat divergence, has 
no interpretation out of equilibrium. 

The order parameter in the present case is a non- 
conserved complex scalar field. Its dynamics, however, 
is not expected to be autonomous even close to the crit- 
ical point as density fluctuations are needed to create 
the ordered phase. Although it has been shown that the 
transition dynamics is mediated by waves [l3 |. momen- 
tum density decays fast due to friction. Therefore, the 
most appropriate description in the theory of dynamical 
critical phenomena is model C, in which a non-conserved 
order parameter is coupled to a conserved non-critical 
density 16[. In this case lid [21 1 and in extensions to 
non-equilibrium dynamics |22| the dynamical exponent 
is predicted to be z = 2 + a/ v^, consistent with the 
measurements if a = 0. 

Conclusions. We have demonstrated that the non- 
equilibrium solid-liquid transition that occurs in a shal- 
low, quasi-two-dimensional granular system can be of ei- 
ther first or second order type depending on the vertical 
height and filling density. This seems counterintuitive, 
because it is widely believed that a solid-liquid phase 
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transition can only be of first order. However, motivated 
by observations inside carbon nanotubes recent molec- 
ular dynamic simulations show that in confined water 
nanofilms the transition to a solid phase can be either of 
first or second order, depending on the filling density (23j . 
In our experiments, for both cases density fluctuations 
do not show strong variations at the transition. On the 
contrary, local order varies strongly, either abruptly in 
the first order type transition, or continuously in the sec- 
ond order type configuration. The continuous transition 
presents critical-like behavior, with exponents consistent 
with model C of dynamical critical phenomena. 
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